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ABSTRACT 

The thermal instability of a couple-stress fluid acted upon by 
uniform vertical rotation and heated from below is investigated. 
Following the linearized stability theory and normal mode 
analysis, the paper through mathematical analysis of the governing 
equations of couple-stress fluid convection with a uniform vertical 
rotation for the case of rigid boundaries shows that the complex 
growth rate £7 of oscillatory perturbations, neutral or unstable for 
all wave numbers, must lie inside a semi-circle 

|oj —T A — (2F;T 6 ), in the right half of a complex <7 -plane, 

where F is the couple-stress parameter, which prescribes the upper 
limits to the complex growth rate of arbitrary oscillatory motions 
of growing amplitude in a rotatory couple-stress fluid heated from 
below. 
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1. INTRODUCTION 

Right from the conceptualizations of 
turbulence, instability of fluid flows is being 
regarded at its root. The thermal instability 
of a fluid layer with maintained adverse 
temperature gradient by heating the 
underside plays an important role in 
Geophysics, interiors of the Earth, 
Oceanography and Atmospheric Physics etc. 
A detailed account of the theoretical and 



experimental study of the onset of Benard 
Convection in Newtonian fluids, under 
varying assumptions of hydrodynamics and 
hydromagnetics, has been given by 
Chandrasekhar 4 . The use of Boussinesq 
approximation has been made throughout, 
which states that the density changes are 
disregarded in all other terms in the equation 
of motion except the external force term. 
Sharma et a/. 8 has considered the effect of 
suspended particles on the onset of Benard 
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convection in hydromagnetics. The fluid has 
been considered to be Newtonian in all 
above studies. With the growing importance 
of non-Newtonian fluids in modern 
technology and industries, the investigations 
on such fluids are desirable. Stokes 12 
proposed and postulated the theory of 
couple-stress fluid. One of the applications 
of couple-stress fluid is its use to the study 
of the mechanism of lubrication of synovial 
joints, which has become the object of 
scientific research. According to the theory 
of Stokes 12 , couple-stresses are found to 
appear in noticeable magnitude in fluids 
having very large molecules. Since the long 
chain hylauronic acid molecules are found as 
additives in synovial fluid, Walicki and 
Walicka 13 modeled synovial fluid as couple- 
stress fluid in human joints. Sharma and 
Thakur 9 have studied the thermal convection 
in couple-stress fluid in porous medium in 
hydromagnetics. An electrically conducting 
couple-stress fluid heated from below in 
porous medium in the presence of uniform 
horizontal magnetic field has been studied 
by Sharma and Sharma 10 . Sharma and 
Sharma 1 and Kumar and Kumar 5 have 
studied the effect of dust particles, magnetic 
field and rotation on couple-stress fluid 
heated from below and for the case of 
stationary convection, found that dust 
particles have destabilizing effect on the 
system, where as the rotation is found to 
have stabilizing effect on the system, 
however couple-stress and magnetic field 
are found to have both stabilizing and 
destabilizing effects under certain 
conditions. 

Banerjee et al. 2 gave a new scheme 
for combining the governing equations of 
thermohaline convection which is shown to 



lead to bounds for the complex growth rate 
of the arbitrary oscillatory perturbations, 
neutral or unstable for all combinations of 
dynamically rigid or free boundaries. 
However no such result existed for non- 
Newtonian fluid configurations, in general 
and for couple-stress fluid configurations, in 
particular. Banyal 3 have characterized the 
non-oscillatory motions in couple-stress 
fluid. 

Keeping in mind the importance of 
non-Newtonian fluids, the present paper is 
an attempt to prescribe the upper limits to 
the complex growth rate of arbitrary 
oscillatory motions of growing amplitude, in 
a layer of incompressible couple-stress fluid 
heated from below in the presence of 
uniform vertical rotation opposite to force 
field of gravity, when the bounding surfaces 
of infinite horizontal extension, at the top 
and bottom of the fluid are rigid. 

2. FORMULATION OF THE PROBLEM 
AND PERTURBATION EQUATIONS 

Considered an infinite, horizontal, 
incompressible couple-stress fluid layer, of 
thickness d, heated from below so that, the 
temperature and density at the bottom 

surface z = 0 are T 0 , p 0 respectively and at 

the upper surface z = d are T d , p d and that a 
uniform adverse temperature gradient 



P 



f 


dT 




V 


dz 


J 



is maintained. The fluid is acted 



upon by a uniform vertical rotation 

fl(0,0, Cl). Let/?, p, T and q(u,v,w) 
denote respectively the density, pressure, 
temperature and velocity of the fluid. Then 
the momentum balance, mass balance 
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equations of the couple-stress fluid The basic motionless solution is 
(Stokes 12 ; Chandrasekhar 4 and Scanlon and 
Segel 6 ) are 



dq_ 

dt 



q y\q 



P „2 L2 



P0 

( 



-Vp+ g 1 + 



Sp_ 
P0 



r _j^_V^ q + 2 
P0 



— > — > 
qxQ. 



(1) 



The equation of state 
p = p 0 [l-a{T-T 0 )], 



(2) 



(3) 



Where the suffix zero refer to the values at 

the reference level z = 0. Here g(0,0,-g) is 

acceleration due to gravity. 

Let c v denote the heat capacity of 

the fluid at constant volume. Assuming that 
the particles and the fluid are in thermal 
equilibrium, the equation of heat conduction 
gives 



Po C v 



Or 

dt 



+ 



— + q.V T = qV 2 T, 
dt ) 



q$ T = J& 2 T , 

V J 



(4) 



The kinematic viscosity V , couple- 
stress viscosity// , thermal diffusivity 



K = 



, and coefficient of thermal 



expansion a are all assumed to be constants. 



9 = (0,0,0) , T = T Q -/3z, n = (o,o,fl) 

and p = p Q {l + a/k). (5) 

Assume small perturbations around 
the basic solution and let dp , dp , 6 and 

q(u,v,w) denote respectively the 
perturbations in density, pressure p, 
temperature T and couple-stress fluid 
velocity (0,0,0) . The change in density dp 

caused mainly by the perturbation 0 in 
temperature is given by 



dp = -ap Q 0 . 



(6) 



Then the linearized perturbation 
equations of the couple-stress fluid becomes 



do 



V--^-V 2 V 2 q+lqxCl 



P0 



dt 



= J3w + /ft 2 0, 



(7) 

(8) 
(9) 



Within the framework of Boussinesq 
approximation, equations (7) and (8), give 



— V'w-ga 
dt 



dx 2 dy 2 ) 



"3? 






( i A 


dt 






P0 

V J 



dz 



P0 



(10) 

(11) 
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where 



Together with (9), where 

V*- 52 + d * + ** a = kd,a = ^ Pl = V -,F = ^ 2 

~dx 2 dy 2 dz 2 and V K PodV 



„ 4Q 2 d 4 d 
T A = — D = — and D m = dD and 

_dv du v & 

— -— denote the z-component of dropping (©) for convenience. Here 

vorticity. Pi = ~~> is the thermal prandtl number, F is 

3. NORMAL MODE ANALYSIS the couple-stress parameter and T A is the 

Taylor number. Substituting W = W @ , 
Analyzing the disturbances into 2 r— 

normal modes, we assume that the q _ P q an( j Z — Z in 
Perturbation quantities are of the form K d 

equations (13), (14) and (15) and dropping 
[w,0, g\- [w(z),©(z),Z(z)] (©) for convenience, in non-dimensional 

Exp{ik x x + ik y y + nt) (12) form becomes, 

Where k x ,k y are the wave numbers along (D 2 -a 2 |<7+fjD 2 -a 2 j -(o 2 -a 2 j w=-rc?&-t a dz ( (16) 
the x and y-directions respectively 

k = {k x + ^ y 2 ) 2 , is the resultant wave 

number and n is the growth rate which is, in 
general, a complex constant. 
Using (12), equations (9), (10) and (11), on 
using (8), in non-dimensional form, become 



l-F]D 2 -a 2 }\ D 2 -a 2 



Z - DW , (17) 

(d 2 - a 2 - p,a)& = -W , (18) 



D 2 -a 2 



gaBd 4 

Where R = , is the thermal Rayleigh 

a + F\D 2 -a 2 Y-(D 2 -a 2 )\ KV 



number. 

gad 2 a 2 0 I — ' Since both the boundaries rigid and 

~ v yjTAdDZ ar£ mamtame( j at constant temperature, the 

perturbations in the temperature are zero at 
i f(d 2 ■J 2 l]f-D 2 j 2 ) gz= ^ a dw (14) the boundaries. The appropriate boundary 
J d conditions with respect to which equations 

(16), (17) and (18) must be solved are 

D 2_ fl 2 _ ) &= _&l w W = DW = 0, 0 = 0 and Z = 0 at z = 0 and 

) k > Z=L (19) 
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Equations (16)-(18), along with Taking complex conjugate on both sides of 

boundary conditions (19), pose an equation (18), we get 
eigenvalue problem for (7 and we wish to 

Characterize <7, when 0 r > 0 . (d 2 - a - p x CT*)&* = -W* , (21) 
We prove the following theorem: 

Theorem: If R) 0, F>0, T A )0, <7 r >0 

and <7 i ^ 0 then the necessary condition for ( 

the existence of non-trivial solution Jw *®dz = - j"©(z) 2 - a 2 - p x a* p*dz > ( 22 ) 

{W,©,Z) of equations (16), (17) and (18) 0 

together with boundary conditions (19) is Also takin § complex conjugate on both sides 

that|cr| 2 <r,-(2^ 6 F). ofequationd^weget 

r\D 2 -a 2 )\(D 2 -a 2 \- o Z =-D\V~ f (23) 



Proof: Multiplying equation (16) by W 
(the complex conjugate of W) throughout Therefore, using (23), we get 
and integrating 

the resulting equation over the vertical i i i r 2 ] 

range of z, we get ]\fDZd^-lD\fzd^]z\(iy l -« 2 j-ffo 2 -a 2 ] -a* W dz , 



o o 



o\W*{d 2 -a 2 )wcIz+F\W*(d 2 -« 2 f Wife- j\V* 
0 0 0 



(24) 



Substituting (22) and (24) in the right hand 

V -a 2 )' Wd z = -Ra 2 )w* &dz -T A )w*DZ dz Side ° f eC l Uati0n ( 20 )' We § et 



1 ( \ 1 ( \^ 1 ( \ 2 

a\W*[D 2 - a 2 \Wdz +F \W *\d 2 - a 2 Wdz - \W *[d 2 -a 2 Wdz 



■■ Ra* 



J©(^D 2 -a 2 - Pl a* y*dz-T A \zljD 2 - a 2 ^- f[d 2 - a 2 ^ - a*\z* dz , 



(25) 



Integrating the terms on both sides of boundary conditions (19), along with (17), 
equation (25) for an appropriate number of we get 
times by making use of the appropriate 
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1 ' 



D 2 W 



* I 1 2 



2 +2a 2 |w| 2 +a 4 |w| 2 |dz =«fl 2 }||D@| 2 +a 2 |©| 2 + pi<7*|0| 
|2 , 2\ 7 \2~ 



dz -T A \ \\DZ\ + a 2 \z\ Idz 

0 ^ ' 

-T A F\\ D 2 Z 2 + 2a 2 \DZ\ 2 +a 4 \z\ 2 ldz-T A (T* \\z\ 2 dz 



Ol 



0 



(26) 



And equating imaginary parts on both sides 
of equation (26), and cancelling <7 ; (^ 0) 
throughout from imaginary part, we get 



j||DH'| 2 + o 2 |iv| 2 lfe + «ii 2 p 1 j|0| 2 A = r A J|z| 2 & • 



(27) 



We first note that since W and Z satisfy 
W(0) = 0 = W(\) and Z(0) = 0 = Z(l) in 
addition to satisfying to governing equations 
and hence we have from the Rayleigh-Ritz 
inequality 7 



\\DW\ 2 dz>x 2 \\w\ 2 dz, 



(28) 



and 



\\DZ\ 2 dz>x 2 \\z\ 2 dz, 



(29) 



Further, for W(0) = 0 = W(l) and 

Z(0) = 0 = Z(l) , Banerjee et al. 1 have 

shown that 

i i 

J]D 2 w| 2 Jz > n 1 j|£W| 2 dz and 



§D 2 z\ 2 dz>x 2 \\DZ\ 2 dz, (30) 



Further, multiplying equation (17) and its 
complex conjugate (23), and integrating by 
parts each term on both sides of the resulting 
equation for an appropriate number of times 
and making use of boundary condition on Z 
namely Z(0) = 0 = Z(l) along with (17), 
we get 



l{o 2 z\ 2 +2a 2 \DZ\ 2 +a 4 \z\ 2 \z + F 2 jj \d 4 z\ 2 + 4a 2 \D 3 z\ 2 + 6a 4 \\d 2 z"\ +4a 6 j\DZ\ 2 +a s \z\ 2 \dz 

0 0 [ 0 0 J 



2F JF 



+ 3a' 



D 2 Z 



■ 3a 4 \\DZ\ 2 + a 6 |z| 2 \dz + 2a r \{dz\ 2 +a 2 \z\ 2 \iz 
o Jo 
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+ 2(7 



,f)\d' 



+ 2a 



D 2 Z 



+ a 4 \z\ 2 \tz + \<r\ 2 J|z| 2 <fe = \\DW\ 2 dz, 



(31) 



Further, by utilizing boundary conditions (19) and equation (17), it follows that 



1 

- \DZ*D 3 Zdz 


< 


i 

\DZ*D 3 Zdz 


0 




0 



i 

\\DZ\ 2 dz 


2 


~1 

\\D 3 z\ 2 dz 


_o 




_o 



(Utilizing Cauchy- Schwartz -inequality), 



71 



1 

if 1 



J]d 2 z| J]d 3 z| & 

.0 J L o 

(Utilizing inequality (30)), 



So that we have 
i 



1 


3 <I 




2 


_0 




_o 





1 1 



2F j]z> 3 z|<fe + |cr| 2 j\z\ 2 dz(j\DW \ 2 dz > ( 34 ) 



0 0 



And on utilizing the inequalities (29), (30) 

Which yields and (33), inequality (34) gives 

11 1 1 1 

§D 3 z\ 2 dz>X 2 ttD 2 z\ 2 dz, (32) \\z\ 2 dz{-t— e —^\\DW\ 2 dz, (35) 

00 o F + \<7\ jo 

Using inequality (29) and (30), inequality Now R > 0 and T A ) 0, utilizing the 

(32) becomes inequalities (35), the equation (27) gives, 

1 i 

^D 3 z\ 2 dz>x 6 \\z\ 2 dz, (33) 



2* b F + \ a y 



\ \DW I dz + a ■ 

0 



<0, (36) 



Now F ) 0 and <T r > 0 , therefore the J \w | 2 d Z + Ra 2 P1 f |© | 2 & 
equation (31) gives, 
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and therefore , we must have 
\a\\T A -(2x 6 F). (37) 
Hence, if 

a r > 0 and O i ± 0 , then 
\ct\ 2 (T a -(2x 6 f). 

And this completes the proof of the theorem. 

4. CONCLUSIONS 

The essential content of the theorem, 
from the point of view of linear stability 
theory is that for the configuration of 
couple-stress fluid of infinite horizontal 
extension heated form below, having top and 
bottom bounding surfaces rigid, in the 
presence of uniform vertical rotation parallel 
to the force field of gravity, the complex 
growth rate of an arbitrary oscillatory 
motions of growing amplitude, lies inside a 

semi-circle in the right half of the a r <J i - 

plane whose centre is at the origin and radius 

is jr A -[2F^ . 

Further, it follows from inequality 
(37) that a sufficient condition for the 
validity of the 'principle of exchange of 
stabilities' in rotatory couple-stress fluid 
convection is that the non-dimensional 
number . Ta <i . It is therefore clear that the 

existence of oscillatory motions of growing 
amplitude in this problem depends crucially 
upon the magnitude of the non-dimensional 



number Ta , in the sense so long as 
o( . Ta <i , no such motions are possible. 
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